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fvq , Abstract. The three-dimensional shapes of thin laminae in such as leaves, wings, flowers 

etc. are driven by the differential growth of these objects. The growth takes place through 
variations in the metric as a function of location in the central plane of the object, and 
also across the thickness of the sheet. In this paper, we provide a rigorous derivation 
of asymptotic theories for the shape of residually strained thin lamina with nontrivial 
curvature, in both the strong and weak curvature regimes, extending hence the discussion 
in |15) . In particular, we derive the precise conditions under which the model equations 
proposed for the description of the deployment of petals during the blooming of a flower 
[21] are valid, and we generalize the classical Foppl-von Karman energy for prestrained 
■ shells with various relative orders of curvature. For weakly curved shells, we show that 

the derived models depend on the specific choice of scaling laws for the shallowness of the 
mid-surface and for the growth tensor with respect to the thickness, and lead to different 
asymptotic theories based on these assumptions. Along the way, we also propose a new 
matching property for second order infinitesimal isometries on convex weakly shallow shells. 
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1. Introduction 



Of late, there has been increasing interest in characterizing the growth-induced morpho- 
genesis of low- dimensional structures such as filaments, laminae and their assemblies, which 
arise routinely in biological systems and their artificial mimics [131 HH El H] • The physical 
basis for morphogenesis is now classical and elegantly presented in the Thompson's opus 
"On growth and form" in terms of a simple principle: differential growth in a body leads 
to residual strains that will, generically, result in changes in shape of the body. Eventually, 
^ ; the growth patterns will themselves be regulated by the strains, so that this principle be 
the basis for the physical self-organization of biological tissues. While such questions lay 
at the interface of biology, chemistry and physics, they are also fundamentally of geometric 
nature. Indeed, they may be characterized in terms of a variation on a classical theme in 
differential geometry - that of embedding a shape with a given metric in a space of possibly 
different dimension [231 [251 [9]. However, the goal now is not only to state the conditions 
when it might be done (or not), but also to constructively postulate the resulting shapes 
in terms of an appropriate mechanical theory. 

The combination of the separation of scales that arises naturally in slender structures 
and the constraints associated with the prescription of growth laws that are functions of 
space (and time), leads to the expectation that the resulting theories ought to be variants 
of classical elastic plate and shell theories such as the Foppl-von Karman or the Donnell- 
Mushtari-Vlasov pp. That this is the case, has been shown for bodies that are initially flat 
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and thin i.e. elastic plates with no initial curvature, using analogies to thermoelasticity 
[22. 20J, perturbation analysis [4, 5J, and rigorous asymptotic analysis [15J. However, most 
laminae are naturally curved in their strain-free configurations, which often is a consequence 
of slow relaxation, perhaps following a previous growth history. Since even infinitesimal 
deformations of a curved shell will potentially violate isometry relative to its rest state, one 
expects that differential growth of such an object will likely lead to a variety of possible 
low dimensional theories depending on the relative size of the metric changes imposed on 
the system. This potential multiplicity of asymptotic theories is of course presaged by a 
similar state of affairs for the derivation of a nonlinear theory of elastic shells [HI IB] • 

Following the study of the nonlinear elasticity for thin plates and shells using a rigorous 
asymptotic theory [TJ EJ [TBI El [TS] and a linearized theory [TH] for residually strained 
Kirchhoff plates [12], we recently provided a rigorous derivation for a nonlinear Foppl-von 
Karman like theory for the growth of naturally flat plates [15] . The resulting theory was 
consistent with that postulated in [20] using a thermoelastic analogy, and derived formally 
from 3d elasticity using perturbation theory [U [5]. Here, we continue and extend the 
discussion in [20j (2TJ [15] concerning the rigorous derivation of asymptotic theories for the 
shape of residually strained thin lamina with nontrivial curvatures. 

As our starting point for a similar theory for growing curved shells, we use the observa- 
tion that it is possible to change the shape of a lamina (a blooming lily petal) by driving it 
via excess growth of the margins relative to the center, rather than via the midrib deforma- 
tion, as conjectured in previous studies [28]. Again, a thermoelastic analogy [22] suggests 
a natural generalization of the Donnell-Mushtari-Vlasov shell theory [1] to growing shells 
[2T] proposed as a mathematical model for blooming, activated from the initial (transverse) 
out-of-plane displacement v o of a petal's mid-surface, with respect to a flat reference config- 
uration. When f o = the equations ( 13. ip reduce to the prestrained von Karman equations 
(12. 7p proposed in [20] and rigorously derived in [T5] from the so-called non-Euclidean elas- 
ticity, where the imposed 3d prestrain is given via a Riemannian metric, whose components 
display the appropriate target (linear) stretching tensor e g (of order 2 in shell's thickness 
h), and the bending tensor K g (of order 1 in h, see ( 12. ip ). 

In this paper we perform a rigorous analysis of the 3d non-Euclidean elasticity model 
and derive reduced-order theories for elastic shells that are weakly strained by differential 
growth. In Section [2] we show that for non-flat mid-surface (of arbitrary curvature and when 
v o 7^ is arbitrarily large), the variationally correct limiting theory ( 12. 3 p is an extension of 
the classical von Karman energy to shells in [16]. In the special case Vq = 0, this energy 
still reduces to the functional whose Euler-Lagrange equations give precisely (12. 7p . 

In Section [31 we show that the equations ( 13. ip can be formally derived from ( 12 .7p by 
pulling back the prestrain tensors e g and K g from the shallow shells whose reference mid- 
surface is given by the out-of-plane displacement hvo, hence comparable with the shell's 
thickness h. We then show (see Theorem 13 .3 j) that ( 13. ip is the Euler-Lagrange equations 
of the variational limit for the 3d nonlinear elastic energies on weakly prestrained shallow 
shells, where the order of the prestrain tensor is appropriately compatible with the order 
of the mid- surface curvature. 

In Section HJ we consider a more general situation, where the referential mid-surface is 
given by the out-of-plane displacement h a vo, for an arbitrary exponent a > 0. When a < 1, 
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we prove in Theorem 14. 1 1 that in the limit the shallowness takes over the prestrain, i.e. there 
is no effect of the tensors e g , K g on the limiting variational lower bound of the energy. This 
energy, given in ( 14. 3 p consists of the linearized Kirchhoff-type [8] functional which minimizes 
the relative bending under the constraint ( 14. 4p that the limiting displacement v and the 
original profile Vq have the same determinants of their Hessians. This amounts to saying 
that the image surfaces of id + fiv^e^ and id + hve^ have the same Gaussian curvatures, up 
to the highest order terms. 

We conjecture that the same energy is also the variational upper bound, when a < 1. 
As a first step towards establishing this natural claim, and in order to shed light on the 
mentioned admissibility constraint, we prove Theorem 14.31 corresponding to a matching 
property which states that v is admissible if and only if it can be matched by equibounded 
higher order displacements to be isometrically equivalent with vq. For results similar in spirit 
but in different contexts see [H [27J, EH HO]- The two essential assumptions of Theorem 14.31 
are: regularity v G C 2 ' 13 , and convexity det V 2 t>o > 0. 

We omit most of the proofs in the paper, since they consist of tedious yet only minor 
from the analytical viewpoint, necessary modifications of the arguments in [TBI [15]. The 
result in Theorem 14.31 is however completely new, hence we present its proof entirely, in 
Section 



2. The prestrained von Karman energy on shells of arbitrary curvature 

Let S be a 2-dimensional surface embedded in M. 3 . We assume that S is compact, con- 
nected, oriented, and of class C 1 ' 1 , and that its boundary dS is the union of finitely many 
(possibly none) Lipschitz continuous curves. Consider a family {S h }h>o of thin shells of 
thickness h around S: 

S h = = x + tn(x); xG-S, -^<t<~j, < /i < /i < 1 

where we use the following notation: n(x) for the unit normal, T X S for the tangent space, 
and n(x) = Vn(x) for the shape operator on S, at a given x G S. The projection onto S 
along n is denoted by tt, so that tt(z) = x for all z = x + tn(x) G S h , and we assume that 
/?, <§C 1 is small enough to have 7r well defined on each S h . 

The instantaneous growth of S h is described by two smooth tensors: e g , K g : S — > M 3x3 : 

(2.1) a h = [a h i3 ] : ^ — -> M 3x3 , a h {x + tn) = Id + h 2 e g {x) + htK g (x). 

We note that the growth tensor a h is as in [201 EE] , but now in a general non-flat geometry 
setting. 

For an elastic body, the energy density W : M 3x3 — > R + , in addition to being C 2 regular 
in a neighborhood of SO (3), is only assumed to satisfy normalization, frame indifference 
and nondegeneracy conditions: 

3c> VF G M 3x3 Vi? G SO(3) W(R) = 0, W(RF) = W(F), 

W{F) > cdist 2 (F,50(3)). 
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where F = Vu is the deformation gradient associated with a mapping u. Consequently, 
the thickness averaged elastic energy of a deformation u h is given by: 

(2.2) l h (u h ) = \( W{Vu h {a h )- 1 ) dz, Vw G W 1,2 (S h , R 3 ). 



(2.3) 



Introducing the asymptotic limit (as h — > 0) of the scaled energies J h then leads, following 
Theorems 12.11 and 12.21 to the variationally correct model for weakly prestrained shells with 
arbitrary large curvature, and corresponds to a nonlinear energy functional X 4 acting on 
admissible limiting pairs (VgV,5gB): 

WgV WBeB l4(V,B)=^J^Q 2 ^x,B-^(A% an -(sjme g ) tan ^ 

+ 2~4 / 22 ( V ( A ") ~ ^ n kn - (sym /tjton). 

Here, the space V consists of first-order infinitesimal isometries on 5* defined by: 
(2.4) V = {1/ G H/ 2,2 (5, M 3 ); r • <9 r V(x) = Va.e. x G 5 Vr G T^} , 

that is those W 2,2 regular displacements for whom the change of metric on S due to the 
deformation id+eV^ is of order e 2 , as e — > 0. Furthermore, for a matrix field A G L 2 (S, M 3x3 ), 
by A tan ( x ) we denote the tangential minor of A at x G S", that is [(A(a;)r)7y] TflG T a: 5- We 
note that the skew-symmetric gradient of V as in ( 12.4ft uniquely determines a H^ 1,2 matrix 
field A : 5 1 — )■ so(3) so that: d T V{x) = A{x)r for all r G T^. Hence, we may alternatively 
write: 

V = [V G W/ 2 ' 2 (S, M 3 ); G H/ ll2 (5, M 3x3 ) Va.e. s G 5 Vr e 1,5 

d T V(x) = A(x)t and A(a;) T = -A(x)}. 

For a plate, that is when S C M 2 , an equivalent analytic characterization for V = 
(V 1 , V 2 , V 3 ) G V is given by: (V l ,V 2 ) = (—u)y,ujx) + (61,62), while the out-of-plane dis- 
placement V 3 G W 2,2 (S, M.) remains unconstrained. 

Moreover, the space B in (I2.3P consists of finite strains: 



B= |L 2 -limsymVw 6 ; w e G W X ' 2 {S, M 3 )|, 



which are the limits of symmetrized gradients of all sequences of displacements on S. Above, 
by symVto(i) we mean a bilinear form on T X S given by: (symVu> (x)t)t] = ^[(d T w(x))r) + 
(d„w(x))r], for all r, rj G T X S. 

It follows (via Korn's inequality) that for a naturally flat plate B = {symVw; w G 
W 1,2 (S, M 2 )}, i.e. B consists precisely of symmetrized gradients of all the in-plane displace- 
ments of S. When S is strictly convex, rotationally symmetric, or developable without flat 
regions, it has been proven in [16] that B = L 2 (S', M 2 ^), i.e. it contains all symmetric 
matrix fields on S with square integrable entries. 

Finally, in (12.31) . the quadratic forms : 



Q 3 (F) = D 2 W(ld)(F,F), Q 2 (x,F tan ) = mm{Q 3 (F); F G M 3x3 , (F - F) tan = 0}. 
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where the form Q3 is defined for F G IR 3x3 , while Q.2(%, •), for a given x G S is defined on 
tangential minors F tan of such matrices. Both forms Q3 and all Q,i{x, •) are positive definite 
and depend only on the symmetric parts of their arguments. 

We now have the following results, stating in particular that the functional X 4 is the 
T-limit [3] of the scaled energies h~ A I h : 

Theorem 2.1. Let a sequence of deformations u h G W 1,2 (S h , M 3 ) satisfy I h {u h ) < Ch A . 
Then there exists proper rotations R h G 5*0(3) and translations c h G IR 3 such that for the 
renormalized deformations: 

y h (x + tn(x)) = (R h ) T u h (x + t^-n) - c h : S h ° — ► IR 3 

ho 

defined on the common thin shell S h ° , the following holds. 

(i) y h converge in W 1 ' 2 (S ho ,M?) to n. 

(ii) The scaled displacements: 

nho/2 

(2.5) V h (x) = h- 1 f y h (x + tn)-x dt 

J -ho/2 

converge (up to a subsequence) in W 1,2 (S, R 3 ) to some V G V. 

(iii) The scaled averaged strains: 

(2.6) B h (x) = h- 1 sjmVV h (x) 

converge (up to a subsequence) weakly in L 2 (S,M. 2x2 ) to a limit B G B. 

(iv) The lower bound holds: 

liminf h- 4 I h {u h ) >1±{V,B). 

Theorem 2.2. For every couple V G V and B G B, there exists a sequence of deformations 
u h e w 1 ' 2 (S h ,R 3 ) such that the following holds: 

(i) The rescaled sequence y h (x + tn) = u h (x + tj^n) converges in W 1,2 (S h °, 1R 3 ) to n. 

(ii) The displacements V as in $2.5\) converge in W l,2 (S, M 3 ) to V. 

(iii) The strains B h as in A2.6\) converge in W 1,2 (S, IR 2x2 ) to B. 

(iv) One has: 

lim h- 4 I h {u h ) =Zt(V,B). 

The proofs follow through a combination of the arguments in [T5] and [TB] , and we leave 
them to the interested reader. We now comment on the functional (12. 3 p and its relation 
with the prestrained von Karman equations for plates. 

Here, in analogy with the theory for naturally flat plates with incompatible strains [To] , 
in (12.11) . we have assumed that the target metric is 2nd order in thickness for the in- 
plane stretching: (sym e g ), and 1st order in the thickness for bending: (sym K g ). Due to 
this particular choice of scalings the limit energy X4 is composed of exactly two terms, 
corresponding to stretching and bending. The argument of the integrand in the first term: 
B — \{A 2 ) tan — (sym e g ) tan is the difference of the second order stretching induced by 
the deformation id + hV + h 2 w h from the target stretching (sym e g ), where V G V and 
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symVio — > B. The argument of the integrand in the second term: (V(An) — AU) ta n ~ 
(sym K g )tan is the difference of the first order bending induced by the same deformation 
from the target bending (sym K g ). 

In general, the second order displacement w can be very oscillatory. Due to the non- 
trivial geometry of the mid-surface, the finite strain space B is usually large. In other 
words, a bound on the L 2 norm of the symmetric gradients symVw h implies only a very 
weak bound on w h , which means that the limiting B can be written only as the symmetric 
gradient of a very weakly regular distribution (not a classical higher order displacement). 

Remark 2.3. In the case of a plate, i.e. when S C IR 2 , the term B — \{A 2 ) tan is reduced to 
the classical case: \ (Vw + (Vu>) T + <g> Vv) , where w and v = V 3 are respectively the 
in-plane and the out-of-plane displacements. The term (V(An) — AI\) tan reduces also to: 
—V 2 v. As shown in [15], under the assumption of W being isotropic, the Euler-Lagrange 
equations of X4 can be then written in terms of the displacement v and the Airy stress 
potential $: 



where S is the Young modulus, B the bending stiffness, v the Poisson ratio (given in terms 
of the Lame constants \i and A), and : 



Equations (12. 1\ are based on a thermoelastic analogy to growth [221 ED] and also derived 
using a formal perturbation theory [1] . 

Remark 2.4. When the mid-surface 5* is elliptic, for any first order isometry V G V one 
can take B e B = L 2 (S, M. 2 *^) such that B — \(A 2 ) tan — (sym e g ) ta n = 0. This implies that 
for any first order isometry, there exists a higher order modification w h for which in the 
limit, the second order target stretching is realized. Thus, the energy reduces to: 



i.e. the bending term to be minimized over the space V. Note that this problem is convex 
(minimizing a convex integral over a linear space V), and hence it admits only one solution 
(up to rigid motions). Following the analysis in [17], we see that for elliptic surfaces, 
all limiting theories for h~^I h under the energy scaling (3 > 2, coincide with the linear 
theory X 4 as above, while the sublinear theory, to be used for studying e.g. buckling, is the 
Kirchhoff-like (nonlinear bending) theory corresponding to /3 = 2 and described in [TP]. 



For the case when the mid-surface of the slender structure is a shell with non-zero cur- 
vature, the following system was introduced in [2T], again using a thermoelastic analogy to 



(2.7) 



A 2 $ = -5(det V 2 v + X g ) 
BA 2 v = [v, $] - BVL g , 





3. The prestrained shallow shell model. 
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growth, as a mathematical model of blooming activated by differential lateral growth from 
an initial non-zero transverse displacement field vq'- 

A 2 $ = -5(det V 2 v - det V 2 v + X g ) 



B(A 2 v - A 2 v ) = [v, $] - Bn 6 



(3.1) 

I / ; : \ ~ i' \ - i • — in (it / < v , 

Above, the notation and the growth sources A 5 and Q g are as in ( 12. 8ft . We will now show 
that (13. ip can be derived from the equations ( 12. 7p . 

For a given out-of-plate displacement v G W 2,2 (fl,M.) PI C 1 (0) on an open, bounded 
subset Q C M. 2 , consider a sequence of surfaces: 

S h = 0h(O), where = (x, /iu (a;)) Vx = (xi,x 2 ) G fi. 

The unit normal vector to Sh at <ph{x) is given by: 

W^x) x 9i9fcWl V 1 + ^ |Vf | 2 
We now consider thin plates = x (— h/2,h/2) and thin shallow shells (Sh) h '- 

(3.2) (S„) h = {4(x,x 3 ); x G fi, x 3 G (-h/2,h/2)}, 
where the Kirchhoff-Love extension of (ph on Q h is given by the formula: 

(3.3) <j> h (x,x 3 ) = 4> h {x) + x 3 n h (x) V(x,x 3 ) G 

Let now e g , K g : — > IR 3 ,*^ be two given smooth, symmetric tensors. Then, 

Proposition 3.1. The system A3. Jj) can be derived from the equations A2. 7| ) by 'pulling back 
the prestrain tensors e g and K g from a sequence of shallow shells {Sh) h generated by the 
vanishing out- of -plane displacements Hvq. 

For each small h, define the growth tensors on (Sh) h by: 

(3.4) q h (Mx) + x 3 n h (x)) = ld + h 2 e g (x) + hx 3 K g (x) \/(x,x 3 )ett h . 
The corresponding metric p h = (q h ) 2 on (Sh) h then satisfies: 

p h {(j) h {x) + x 3 n h (x)) = Id + 2h 2 e g (x) + 2hx 3 K g {x) + 0(h 3 ). 
Lemma 3.2. The pull-back of the metric p h through (ph satisfies: 
V(x,x 3 ) en h g h (x,x 3 ) = (V4> h ) T (p h °4>h)(V4> h ) 

= Id + h 2 (2e g (x) + (Vv (a;) ® Vvo(x))*) + 2hx 3 (k 9 (x) - (V 2 v {x))*) + 0(h 3 ), 



where F* G M 3x3 denotes the matrix whose only non-zero entries are in its 2 x 2 minor 
given by F G M 2x2 . 

Proof. By a direct calculation, we obtain: 

dx4> h = (l - x^d^vo, -x 3 hd 2 2 v , hdxv ) + 0{h 3 ), 
d 2 (t>h = ( - x 3 hdf 2 v , 1 - x 3 hd 2 2 v , hd 2 v ) + 0(h 3 ), 

d 3 k = n h ={- hd lVo , -hd 2 v , 1 - ^ 2 |V^D + 0(h 3 ). 
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Hence: 

(V^) T (V^) = Id - 2x 3 h(V 2 v )* + h 2 (Vv ® Vv )* + 0(h 3 ) 
{V4> h ) T (2h 2 e 9 + 2hx 3 K g )(V$ h ) = 2h 2 t g + 2hx 3 n g + 0(h 3 ), 
in view of V<ph = Id + 0(h), and the result follows. 
Proof of Proposition 13.11 

By Lemma 13.21 we see that the growth tensor on Q h has the form: 



(3.5) a h = ^ = ld + h 2 (e g +^(Vv Q ®Vv )^ + hx 3 ^K g - (V 2 v )^ + 0(h 3 ). 

Applying (12. 8p to the modified strain and curvature in a h , to the leading order, we obtain: 
X g (v ) = curFcurl ((e 5 ) 2X 2 + ^Vw ® Vw ) = A 9 + det V 2 v 
£l g (v Q ) = div T div ((k b - V 2 v ) 2 x2 + v cof (k 9 - V 2 v ) 2X 2^ 



where the last equality follows from div cof V 2 t> o = 0. Consequently, (12 .7p for the growth 
tensor (13. 5 j) becomes exactly (13. ip . ■ 

We will now show that (13.11) can be rigorously derived as the Euler- Lagrange equations 
for the T-limit functional of the scaled 3d nonlinear elastic energies on weakly prestrained 
shallow shells. Let v h = u h o <p h g W 1,2 (£l , M 3 ), by means of diffeomorphisms <ph as in 
(13.31) . By a simple change of variables, we see that: 

I h (u h ) = \f^w({Vv h ){Vhr\<l h °hT 1 ) -detV^ dx 

= \ [ wtiVv* 1 )^)- 1 ) ■ detV4> h dx, 
h Jn h v J 

where b h = (q h o 0^) V<^ satisfies 

(b h ) T b h = g h 

and therefore by the polar decomposition of matrices 

b h = R(x,x 3 )a h ontt h 

for some R(x, x 3 ) G SO (3) and the symmetric growth tensor a h is given by (13. 5p . 
For W isotropic, it follows that: 



(3.6) 



I h (u h ) = \ f w((Vv h )(a h )- 1 R(x)- 1 ) -detV4> h dx 
h Jnh \ J 

= ~J w((Vv h )(a h )-^ ■ (1 + 0(h)) dx. 



In other words, heuristically, and modulo the change of variable <ph, the problem reduces 
to the study of deformations of the flat domain Q h with the prestrain tensor a h . Indeed, 
by exactly the same analysis as in [15], Theorems 1.2 and 1.3, we obtain: 



(3. 
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Theorem 3.3. Assume that u h G W 1 ' 2 ((S h ) h ,R 3 ) satisfies I h (u h ) < Ch 4 . Then there 
exists proper rotations R h G 50(3) and translations c h G R 3 such that for the normalized 
deformations: 

(3.7) y h (x,t) = {R h ) T {u h o<t) h )(x,ht)-c h : Q 1 — > R 3 

defined by means of $3.3\) on the common domain Q 1 = Q x (—1/2, 1/2) the following holds: 

(i) y h (x,t) converge in W 1 ' 2 ^ 1 ,M. 3 ) to x. 

(ii) the scaled displacements V h (x) = h~ 4 j-^y 2 y h (x,t) — x dt converge (up to a subse- 
quence) in W l ' 2 (Q,R 3 ) to the vector field of the form (0,0, v) T and v G W 2,2 (Q, R). 

(iii) The scaled in-plane displacements h~ x V^ an converge (up to a subsequence) weakly in 
W 1 ' 2 to we W 1 ' 2 ^^ 2 ). 

(iv) Moreover: liminf^o h~ 4 I h {u h ) >T g ^ Vo (w,v) where: 

T g>vo (w, v) = ~ J Q.2 ^symVw + ^Vv ® V« - ^Vw ® Vw - ( £3)2x2^ 

+ 2l / 22 ~ V2 ^ + ^ 2x2 ) • 

The Euler-Lagrange equations ofZ g ^ are precisely ( fff. 

In the same manner, applying the proof of Theorem 1.4 of [15] to ( I3.6p . yields: 

Theorem 3.4. For every v G W 2,2 (Q,R) and w G W 1,2 (f],R 2 ), there exists a sequence of 
deformations u h G W l,2 ((Sh) h , R 3 ) suc/i £/ia£: 

(i) T/ie sequence y h (x, t) = u h (x + hvo(x)e 3 + htn h (x)) converges in W 1,2 (Q l , R 3 ) to x. 

(ii) The displacements V h as in (ii) Theorem \3. 8\ converge in W 1,2 to (0,0, v). 

(iii) The in-plane displacements /i _1 V^ n converge in W 1 ' 2 to w. 

(iv) \im h ^ h- 4 I h (u h ) =I 9tV0 (w,v). 

Remark 3.5. Comparing functional (13. 8 p with ( 12. 3D . note that the space Vi(Sh) of first- 
order infinitesimal isometries on Sh is made of displacements V : Sh — >■ R 3 of the form: 

V(0 fc (a;)) = (-/ifi(a;),-/if2(aj),f3) Vs G fi, 

such that (^1,^2,^3) G Jy 2 ' 2 (f2, R 3 ) and symV(t>i, t^) = sym(Vt>3 <S> V^o). 
Indeed, the condition symVV^ = on Sh becomes: 

= ^{ d i( v ° <Ph) ■ dj(p h + dj(V o cj) h ) ■ di4> h ) = -h sym[V(vi,v 2 ) - Vt> 3 <g> V«o]y. 
We also see that v 3 can be completed by (vi,v 2 ) to V G Vi(5^) as in (13. 9 P only if: 
(3.10) cofV 2 w : V 2 w 3 = 0, 

the latter being also a sufficient condition when 5 is simply connected. This follows from: 

curl T curl^sym(Vf3 <g> Vt>o)) = curl T curl ( Vt> 3 <S> Vt>oj 

= 5 22 (5if3 ■ div ) + d n (d 2 v 3 ■ d 2 v ) - di 2 (d x v 3 ■ d 2 v + d 2 v 3 ■ d x v Q ) 
= - (dim ■ d 22 v + d 22 v 3 ■ d u v - 2d 12 v 3 ■ d 12 v ) = -cofV 2 w : V 2 w 3 . 
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Hence, we see that the admissible out-of-plane displacements V3 relevant in (I2.3p . must 
obey for the least the constraint (I3.10p . in contrast with the unconstrained 2-scale limiting 
theory (ET5)) . 

4. Shallow shells of lower orders 

More generally, for any exponent a > we may consider the following energies of defor- 
mations on weakly prestrained shallow shells: 



(4.1) l\ u h ) = \ f WdVu^)- 1 ) Vu h G W 1 ' 2 ((S h ) h ,M 3 ), 

h hs h ) h 



with the growth tensor q h given by (|3.4p . on thin shells of the form (13 .2p around the shallow 
mid-surface: 

S h = Mty; <t>h{x) = (x,h a v (x)); v G C 2 >P(Q,R). 

According to our notation, n h (x) stands for the unit normal vector to Sh at <fih{x), i.e.: 
n h (x) = {dx<f> h {x) x d^hix^/ld^hix) x d^^x)]. 

When a > 1, it is clear from the previous calculations that the prestrain takes over the 
effect of shallowness and hence the limiting theory is that derived in [T3], coinciding with 
results of Theorem 13.31 for the case v = and yielding the Euler-Lagrange equations (12. 7p . 

Assume now that < a < 1. By a straightforward calculation, we see that I h (u h ) can 
be written as in (13.61) . with: 

a h = Id + hi 2a (yv ® V^o)* - h a x 3 {V 2 v )* + o{h 2a ) + x 3 o{h a ) 

where the uniform quantities in o(h 2a ), o(h a ) are independent of x%. Following the proof 
of Theorem 1.3 in Ho], we obtain: 



Theorem 4.1. Assume that u h G W l ' 2 ({S h ) h , M 3 ) satisfies I h {u h ) < Ch 2a+2 , where I h is 
given in ( 4-l\ ) and < a < 1. Then there exists R h G 5*0(3) and translations c h G M. 3 such 
that for y h as in \3. 7| ), we have: 

(i) y h (x,t) converge in W 1 ' 2 ^ 1 , M 3 ) to x. 

(ii) the scaled displacements V h (x) = h~ a j-^y 2 y h (x,t) — x dt converge (up to a subse- 
quence) in W l ' 2 {yt^) to (0,0, v) T where v G W 2 ' 2 (tt,R) and: 

(4.2) detV 2 ^ = detVV 

(iii) Moreover: limmfh^oh~( 2a+2 ^I h (u h ) > Z Vo (w,v) where: 

(4-3) l V0 (v) = ±^Q 2 (v 2 v-V 2 v 

The constraint in the assertion (ii) follows by observing that /i^symVV^ converges in 
L 2 (VL) to |(Vu Q ® Vf - Vv ® Vv) which implies (1472]) . 

We now conjecture that the linearized Kirchhoff-like energy (14. 3 p (14.21) is also the exact 
T-limit of the rescaled energies h~^ 2a+2 ^ I h (u h ) when Q is simply connected. If f2 is not 
simply connected, one must replace (14.21) with a more general variant which states that 
Vt>o ® Vt>o — Vf g) Vf is a symmetric gradient. As a first step towards the rigorous proof 
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of this statement, and to shed light on the constraint condition (14.21) . we now derive a 
matching property for weakly shallow shells. Let us first explain the meaning of "matching 
property" and put our analysis in a broader context. 

Remark 4.2. In [18], the authors put forward a conjecture regarding existence of infinitely 
many small slope shell theories (with no prestrain) each valid for a corresponding range of 
energy scalings. This conjecture is based on formal asymptotic expansions and it is in accor- 
dance with the previously obtained results for plates and shells. It predicts the form of the 
2-dimensional limit energy functional, and identifies the space of admissible deformations 
as infinitesimal isometries of a given integer order N > determined by the magnitude of 
the elastic energy. Hence, the influence of shells geometry on its qualitative response to 
an external force, i.e. the shells rigidity, is reflected in a hierarchy of functional spaces of 
isometries (and infinitesimal isometries) arising as constraints of the derived theories. 

In certain cases, a given iVth order infinitesimal isometry can be modified by higher 
order corrections to yield an infinitesimal isometry of order M > N, a property to which 
we refer to by "matching property of infinitesimal isometries". This feature, combined with 
certain density results for spaces of isometries, cause the theories corresponding to orders of 
infinitesimal isometries between iV and M, to collapse all into one and the same theory. The 
examples of such behavior are observed for plates [8], where any second order infinitesimal 
isometry can be matched to an exact isometry (M = oo), for convex shells [17] , where any 
first order infinitesimal isometry satisfies the same property, and for non-flat developable 
surfaces J27J [10] where first order isometries can be matched to higher order isometries. The 
effects of these geometric properties on the elasticity of thin films are drastic. A plate whose 
boundary is at least partially free possesses three types of small-slope theories: the linear 
theory, the von Karman theory and the linearized Kirchhoff theory, whereas the only small 
slope theory for a convex shell with free boundary is the linear theory [17J: a convex shell 
transitions directly from the linear regime to the fully nonlinear bending one if the applied 
forces are adequately increased. In other words, while the von Karman theory describes the 
buckling of thin plates at a body force magnitude of order thickness-cubed, the equivalent, 
variationally correct theory for buckling of elliptic shells is the purely nonlinear bending 
theory which comes only into effect when the body forces reach to a magnitude of order 
thickness-squared. 

In the present context, we have the following: 

Theorem 4.3. Assume that Q is simply connected and let t> G C 2 '^(fl, M.) with det V 2 t>o > 
c> 0. Let v e C 2 ' p {n,R) satisfy: 

(4.4) det V 2 v = det V 2 v in fi. 
Then there exists a sequence Wh € C 2,/3 (S7, IR. 3 ) such that: 

(4.5) V7i > V '(id + hve 3 + h 2 w h ) T V (id + hve 3 + h 2 w h ) = V(id + hv e 3 ) T V(id + hv e 3 ) 
and sup ||w/ l || C 2, / 3 < +oo. 
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Before giving the proof in Section El we present a brief discussion. Writing Wh = Wh,tan + 
w\e% where Wh,tan{ x ) G M? and w\ G R, equation (14.51) becomes: 

Id + /i 2 (2symVu^ to n + Vt> <g> Vu) + 2/i 3 sym(Vt> g> Vw^) 

' U)) + h\{Vw hMn ) T Vw hMn + Vw 3 h ® Vto^) = Id + /i 2 Vf ® Vv - 

Recall that (14.41) is equivalent to curl T curl(Vt> ® Vf — Vt>o ® V^o) = 0, and hence to: 
Vw (g> Vf — Vt>o ® Vt>o = symVwj for some w G C 2 ^(ft,M 2 ) since f2 is simply connected. 
Hence the constraint (I4.4p is necessary and sufficient for matching the lowest order (h 2 ) 
terms in ( 14. 6p . 

Our result states that actually it is possible to perturb w by an equibounded 3d displace- 
ment Wh — w so that the full equality ( 14. 6 p holds. A natural way for proving this is by 
implicit function theorem. Indeed, this is how we proceed, and the ellipticity assumption 
det V 2 f o > is precisely a sufficient condition for the invertibility of the implicit derivative 
£(p) : Cj'^(fi,R) — > C?' P [Q,R), C{p) = -cofV 2 f : V 2 p where p is the variation in w 3 h . An 
extra argument for the uniform boundedness of Wh.tan in C 2 ' 13 concludes the proof. 

Remark 4.4. The condition (14. 5 j) means that each deformation u h : Sh — >■ R 3 of a surface 
Sh = {x + hvo(x)e 3 ; x G Q}, given by u h (x + /iu (x)e 3 ) = x + hv(x)e 3 + h 2 Wh(x) is an 
isometry of Sh- In other words, the pull-back metrics from the Euclidean metric of Sh 
and of u h (Sh) = {x + hv{x)e 3 + h 2 Wh(x); x G Q} coincide. Hence Theorem 14.31 asserts 
that if two convex out-of-plane displacements of first order have the same determinants of 
Hessians, then they can be matched by a family of equibounded higher order displacements 
(the fields Wh) to be isometrically equivalent. When the parameter h is replaced by h a , 
this result can be also viewed in the context of shallow shells. For other results concerning 
matching of isometries see [5J Theorem 7], [T71 Theorem 1.1] and [TUJ Theorem 3.1] (which 
is comparable with [27J Lemma 3.3] and the remark which follows therein). 

5. A MATCHING PROPERTY ON CONVEX SHALLOW SHELLS: PROOF OF THEOREM 14.31 

By a direct calculation, (14.51) is equivalent to: 

(5.1) V(id + h 2 w h ,tan) T V(id + h 2 w h ,tan) = Id + h 2 Vv ® Vv o - h 2 ( Vv + Vz h ) <g> ( Vv + Vz h ) , 

where Whjan C 2,l3 (Q,M. 2 ) and z h = hw\ G C 2,/3 (fi,lR) so that Wh = Wh.tan + w h e 3 * s the 
required correction in (14. 5p . 

1. We shall first find the formula for the Gaussian curvature of the 2d metric in the right 
hand side of (15.11) : 

(5.2) g h (z h ) = Id + h 2 Vv <g> Vv - h 2 (Vv + Vz h ) <8> (Vv + Vz h ). 

Lemma 5.1. Let v ,v G C 2,/3 (fi,lR) and consider the C 1,13 regular metrics on Q of the type: 

g = [gij}ij=i,2 = Id + h 2 (Vv <8> Vv - Vv 1 <g> Vvi). 

Then, for any h > small, the Gaussian curvature n{g) of g is C ^ regular and it is given 
by the formula: 

(5.3) K (g) = h 2 



det(V 2 v 



(l-h^d^d^))' 



(1 - WigvdiVodjVoWdetg (l - ^V^] 2 )' 



■detVV 
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where the Christoffel symbols of g, the inverse of g, and its determinant are: 



1 



(5-4) rj = -g kl (d j9a + d i9jl - d j9ij ) , 

(5.5) = m = ^cof[,,], 

det g = 1 - h^Wvo) 1 ■ Vwi| 2 + h 2 {\Vv \ 2 - \V Vl \ 2 ). 
Proof. Assume first that vq and v\ are in fact smooth. By Lemma 2.1.2 in [9], we have: 

detW 



«(ld - h 2 Vvi ® Vvi) = -/i 

i 

K (#) 



1 -h 2 \v Vl \ 2 ) 2 

h 2 det(V 2 v - [rfA^y 



(l - tfffidiVodjVo)) detg 



(l-h^diVodjVo))' 

Since the two metrics above are equal, the formula ( j5.3{) follows directly. The formula for 
detg is obtained by a direct calculation, via det (A + B) — detA + coL4 : B + det 5, valid 
for 2 x 2 matrices A, B. 

In the general case when vq,Vi are only C 2 ^ regular, one may approximate them by 
smooth sequences Vq,Vi- Then, each K n = k (Id + h 2 ( ®v%- V< g> <)) is given by the 
formula in (I5.3p . and the sequence K n converges in C°'° to the right hand side in (15. 3p . On 
the other hand, K n converges in D'(Q) to n(g), which follows from the definition of Gauss 
curvature k = i? 1212 /detg. Hence the lemma is proven. ■ 

2. Applying Lemma I5TT1 to v\ = v + Zh, we now see that for small h, the Gauss curvature 
of metric gt{zh) vanishes: 

(5.6) K{g h {z h )) = 
if and only if: 

(5.7) $(h,z h ) = 0, 
where: 

$(/*, z)=(l- h 2 \Vv + Vz\ 2 ) 2 det(V 2 v - [T^dkVolij) 



- (1 - tftf'diVodjVo)) d(h, z)det(V 2 v + V 2 z) 

Here: 



d(h,z) = l-h A \(Vv ) x - V(v + z)\ 2 + h 2 (\Vv \ 2 - \Vv + Vz\ 2 ) 
and r^- and g^ are given by ( 15. 4p and ( 15.51) for the metric g = Id + h 2 Vv <8> Vt>o — h 2 (Vv + 
V z) ® (Vf + Vz). We shall consider: 

$ : (-e,e) x C 2 '^(O,R) — > C 0,/J (fi,R) 

and seek for solutions G Cq (f2, M) of (15. 7p with zero boundary data. It is elementary to 
check that $ is continuously Frechet different iable at (0, 0) and that 

$(0, 0) = detV 2 v - detV 2 w = 0. 



is a 
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Moreover, the partial Frechet derivative C = dQ/dz(0,0) : Cq ,i3 (Q,M.) - 
linear continuous operator of the form: 

Wz e C 2 / C(z) = lim -$(0, ez) = lim - ( det V 2 v - det(V 2 v + eV 2 z)) 

e— >o e e 

= lim - ( - e 2 det V 2 z - ecofV 2 v : V 2 z) = -cofV 2 t> : V 2 z. 
e v 7 

Clearly, L is invertible to a continuous linear operator, because of the uniform ellipticity of 
the matrix field V 2 t> which follows from the convexity assumption of det V 2 t> = det V 2 i?o 
being strictly positive. Thus, invoking the implicit function theorem we obtain the solution 
operator: 

such that Zh = Z(h) satisfies ( 15. 7p . Moreover, Z is differentiable at h — and: 



Z'(0) = £- 1 o^(0,0))=0, 



because: 
^(0,0)=(cofV 2 ^o): 

Consequently: 
(5.8) 



d_ 

dh 



J- det[T* 9 fc « ] - - f J^(0, 0)) det V 2 v = 0. 



h 



\zh\\c 2 >ti — > as /i — > 0. 



2 /3 

3. In conclusion, we have so far obtained a uniformly bounded sequence of C out-of- 
plane displacements w\ = z h /h such that the Gauss curvature (15.61) of the metric gh{zh) 
in the right hand side of (15. ip is 0. By the result in [23] it follows that for each small h 
there exists exactly one (up to fixed rotations) orientation preserving isometric immersion 
0, eC 2 (Q,R 2 )oig h (z h ): 

(5.9) V0£V0 h = g h (z h ); det V0 h > 0. 

What remains to be proven is that, in fact, 4>h = id + h 2 Wh,tan with some Wh,tan uniformly 
bounded in C 2 ^(n,R 2 ). 

It is a well known calculation (see [2j [23] ) that (15. 9p implies (is actually equivalent to): 



(5.10) 



v 2 0, - [m 



0. 



where are the Christoffel symbols (15 .4p of the metric g = gh(zh) in (I5.2p . By (I5.9P 
||V0fe||L°° < C, and by (I5.10p ||V 2 0/i||l°o < C, hence ||0/j|| C 2,^ < C. But r| are uniformly 
bounded (with respect to small h) in C 0,13 so by (I5.10P ||V 2 



< C and thus: 



< C. 



Note now that II V 



Hence, by f l5~T0|) : 
(5.11) 



\C0,f> 



< Ch 2 in view of the particular structure of the metrics gh{zh)- 



|VVaIIc°./> < Ch 2 . 
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Therefore, for some A h G M 2x2 we have: 

(5.12) || V0 h - A h \\ C i,e <Ch 2 . 

We now prove that the matrix Ah in the inequality above can be chosen as a rotation and 
hence, without loss of generality, Ah = Id. For each x G Vt there holds: 

(5.13) dist(A h ,SO(3)) < \A h - V0^(x)|+dist(V0 ft (x),SO(3)). 

To evaluate the last term above, write: -\/V0^(x)V0h(x) = QDQ T for some Q G £0(3) 
and .D = diag(Ai, A2) with Ai, A2 > 0. Since det V(ph > 0, it follows by polar decomposition 
theorem that: 

dist(V0 h (x),£O(3)) = \^V4>l(x)V<j> h (x) -W| < C7|D-Id| 

= CmaxjlAi - 1|} < Cmax{|A 2 - 1|} < C|D 2 - Id| 

i i 

= C\Q T V<f h (x)V<p h (x)Q - Id| < C\V<l>lV<i> h {x) - Id| < C7i 2 . 

By the above and (15.131) . ( 15.12ft we see that dist(A/ l , £0(3)) < Ch 2 . Hence, without loss of 
generality, || — Id|| c i„s < Ch 2 and: 

Uh-id\\ c2 ,, <Ch 2 . 

Consequently, (f>h = id + h 2 Wh,tan with Hu^fanllc 2 ./ 3 < C. This concludes the proof of 
Theorem 14.31 in view of (15.1 ft which is equivalent to ( 14.5ft . ■ 

Remark 5.2. In view of the properties of the elliptic Monge- Ampere equation we expect 
that the C 2,13 scalar fields v are dense in the set of the W 2 ' 2 fields with a prescribed, strictly 
positive det V 2 . With such a result at hand, it would indeed follow that for convex shells, 
the linearized Kirchhoff-type energy ( 14. 3 p is the rigorous variational limiting theory in the 
present context of weakly shallow shells, in the same spirit as the corresponding matching 
and density of first order isometries on convex shells [17] resulted in that the only small 
slope theory for an elastic convex shell is the linear theory. The latter problems, when 
posed for surfaces of arbitrary geometry, are more difficult. One could hope to be able to 
prove similar results in case S is a strictly hyperbolic surface. In its generality however 
they reduce to the study of nonlinear PDEs of mixed type for which not so many suitable 
methods are at hand. 

6. Discussion 

Our analysis has shown how one can rigorously derive a general theory of shells with 
residual strain that might arise from relative growth, inhomogeneous swelling, plasticity 
etc. In fact, there are many such theories; each is a consequence of the scalings for the 
magnitude of the shallowness of the shell relative to the magnitude of the strain incompat- 
ibility induced by the metric and curvature growth tensors. Thus, not only do we recover 
a recently postulated model [21], but furthermore we also find a degenerate model in H] 
that corresponds to a constrained shell theory with no contributions from growth. For 
simplicity of presentation, our analysis here is limited to a subset of vast possible scenarios 
since we fix the scalings of the incompatible strain tensors (12. ip or (13. 4p throughout the 
paper. This situation might be contrasted with the case of the equations for the growth 
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of an elastic plate (15] where there seems to be just a single natural choice for the scalings 
which result in a generalized Foppl-von Karman theory, but is similar to the situation for 
the weakly nonlinear theory of elastic shells where there are multiple potential scalings and 
corresponding theories. 

A natural generalization of our results would be to allow for different scaling regimes for 
the growth tensors in search of other possible limiting theories. The same mathematical 
approach can be adapted to these other regimes with variants of Theorem 14.31 as an es- 
sential tool for a rigorous justification of other possible models. Overall, there are three 
independent parameters associated with scaling of the shallowness, and the two incompati- 
ble strains, characterized in terms of their dependence on the thickness in the form h a ; the 
resulting theories depend on the choice of scalings for the three parameters. Thus, there 
is no single correct model in general, but of course when dealing with a concrete situation, 
a choice of particular scaling parameters for the relative magnitude of the thickness, the 
shallowness and the differential growth determine the effective theory. 
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